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Rigorous upper limits on the vertical heat transport in two dimensional Rayleigh-Benard convec- 
tion between stress-free isothermal boundaries are derived from the Boussinesq approximation of 
the Navier-Stokes equations. The Nusselt number Nu is bounded in terms of the Rayleigh number 
Ra according to Nu < 0.2295 Ra^^^^ uniformly in the Prandtl number Pr. This Nusselt number 
scaling challenges some theoretical arguments regarding the asymptotic high Rayleigh number heat 
transport by turbulent convection. 
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Rayleigh-Benard convection is the buoyancy-driven 
flow of a fluid heated from below and cooled from above. 
It is important for a variety of systems in the engi- 
neering, geophysical, and astrophysical sciences, and it 
has long served as a fundamental paradigm of nonlin- 
ear science, chaos, and pattern formation. Indeed, the 
Boussinesq approximation to the Navicr-Stokcs equa- 
tions with the boundary conditions analyzed in this paper 
was Rayleigh's original model for calculating conditions 
for onset [Tj, it is the basis of the Lorenz equations [5], 
and it formed the foundation of developments in the mod- 
ern mathematical theory of amplitude [3] and modulation 
[1] equations. Most recently Rayleigh-Benard convection 
has been the focus of a large body of experimental, com- 
putational, theoretical, and mathematical research aimed 
at characterizing the fully turbulent dynamics for appli- 
cation in geophysical and astrophysical regimes [S]. 

Convective fluid flow increases vertical heat transport 
beyond the purely conductive flux. The dimensionless 
enhancement factor, the Nusselt number Nu, is both of 
fundamental interest for applications and the natural and 
widely recognized measure of the intensity and effective- 
ness of the motion. The most basic question for Rayleigh- 
Benard convection is the dependence of Nu on (i) the 
strength of the thermal forcing, commonly expressed in 
terms of a dimensionless Rayleigh number Ra, (ii) the 
material properties of the fluid, which within the Boussi- 
nesq approximation is set by the dimensionless Prandtl 
number Pr, the ratio of the fluid's momentum and ther- 
mal diffusion coefflcients, (iii) the geometry, typically the 
aspect ratio of the container, and (iv) the boundary con- 
ditions. The connection between these variables is gen- 
erally complex and often not even unique, but in the 
"ultimate" high Rayleigh number regime when the flow 
is turbulent, the presumed functional relation between 
the Nu, Pr, and Ra is Nu - Pr'^ Ra*^. 
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Experiments and simulations with Pr = 0(1) and 
no-slip boundary conditions suggest a scaling exponent 
0-27 < /3 < 0.40 at the highest avaflable Ra [51i]. Vari- 
ous theories suggest (modulo possible logarithmic correc- 
tions) that Nu ~ Pr^/^ Ra^/^ as Ra ^- oo [7HS]. Rigorous 
analyses of the Boussinesq model with no-slip velocity 
and isothermal (flxed temperature) [lOl [11] or flxed heat 
flux [12j or mixed temperature [13j boundary conditions 
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yield upper bounds of the form Nu < cRa ' with pref- 
actors < c < oo independent of Pr, so /? = ^ and 
7 = ^ cannot both hold for very large Pr. The Nu-Ra 
relation is certainly different at Pr — oo where theory 
suggests [Hj and analysis proves [IS] (modulo possible 
logarithmic corrections) that Nu < Ra^^'^. 

Two dimensional Rayleigh-Benard convection displays 
many of the physical and turbulent transport features of 
three dimensional convection and has long been utilized 
as a test-bed for theoretical concepts PJl [13 • The ef- 
fect of free-slip (no-stress) velocity boundary conditions 
on developed turbulent convection has largely been unex- 
plored although we note that the rigorous scaling bound 
reported here was anticipated by recent numerical and 
perturbative investigations of transport limits for finite 
[TH] and infinite i20 Prandtl numbers. This letter 
bridges that gap with a proof that Nu < 0.2295 Ra^/^"^ 
uniformly in < Pr < oo for the Boussinesq model 
in two spatial dimensions with fixed temperature and 
free-slip boundaries. This result refutes predictions of 
a Nu ~ Ra}^'^ ultimate regime insofar as the theoreti- 
cal arguments do not refer specifically to the boundary 
conditions or the spatial dimension. This issue is dis- 
cussed further in the conclusion section at the end of 
the paper. Meanwhile the proof of the bound is pre- 
sented in sufficient detail immediately below for moti- 
vated readers to reproduce the calculation in its entirety. 
The key new idea used to derive the result emerged from 
intuition developed in numerical studies of upper bounds 
[181 [T9] : implement and exploit the bulk averaged en- 
strophy balance available for two-dimensional flows with 
free-slip boundaries to decrease the upper bound. 
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FIG. 1: Geometry for the 2d stress-free convection problem. 
Boundary conditions for T, u, v, and the vorticity cj at the 
isothermal no-slip vertical boundaries are shown. All these 
variables as well as the pressure p are periodic in the horizon- 
tal direction with period F. 



The dimensionless equations of motion for the Boussi- 
nesq approximation are 
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where the Prandtl number Vi ~ v / k in the ratio of the 
fluid's kinematic viscosity v to its thermal diffusivity k, 
and the Rayleigh number Ra = gaAT/i^/i/K where 
g is the acceleration of gravity, a is the fluid's ther- 
mal expansion coefficient, and AT is the imposed tem- 
perature drop across the layer of thickness h. Lengths 
are measured in units of h, time in units of /i^/k, and 
temperature in units of AT. The velocity vector field 
u{x,y,t) = iu{x,y,t)+jv{x,y,t) satisfies no-penetration 
and free-slip (stress- free) boundary conditions, and the 
temperature field T(x,y,t) is isothermal on the vertical 
boundaries at y = and y = 1 as shown in Fig. [T] 
All dependent variables, u, v, T, and the pressure field 
p{x,y,t), are periodic in the horizontal direction x with 
period F (the aspect ratio). 

Taking the curl of ([T]) one obtains the evolution equa- 
tion for the scalar vorticity uj — dv/dx — du/dy, 
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The boundary conditions on u and v imply that lj — Q 
on the vertical boundaries at y = and y = 1. 

The goal of the analysis is to use the equations of mo- 
tion to derive upper bounds on the Nusselt number de- 
fined as Nu = 1 + (fT), where (•) represents the spatial 
and long time average, in terms of Ra, Pr, and F. Toward 
this end we utilize the background method |21| . a math- 
ematical device introduced by Hopf to establish the ex- 
istence of weak solutions to the Navier-Stokes equations 
in bounded domains |22j . For convection problems the 
background method involves decomposing the tempera- 
ture field into a background profile T(y) which satisfies 



the vertical boundary conditions (t(0) = 1 and r(l) = 0) 
and a perturbation term 6{x,y,t) satisfying correspond- 
ing homogeneous boundary conditions {9{x, 0,t) = = 
0ix,l,t)) so that T{x,y,t) ^ T{y) + 9{x,y,t) [TT. Im- 
plementing this decomposition the temperature equation 
([3]) implii 
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Then the equations of motion together with the boundary 
conditions and the background decomposition imply 
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where || • II2 is the L norm on the spatial domain and 
the elementary identity ||Vm||2 — H^Hi was used in 

It is well-known that the equations of motion imply 
Nu = (|Vr|2) [inilll]. Thus, given coefficients a and 
b with precise values to be determined, combining (|6][9]) 
according to 



b 

Ra 



Ra3/2 



^ + 2x^ + ^, (10) 



applying the long time average — remarking that it can be 
shown within the background method that the time aver- 
ages of the time derivatives vanish [Ulini] — and dividing 
by F, the Nusselt number is expressed 
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Hence if we can choose the background profile T(y) and 
coefficients a > and < 6 < 1 so that Q > for all 
relevant 9, lo and w, then the first term on the right hand 
side of ([n]) is an upper bound on Nu. For the problem 
at hand we may use the piece-wise linear profile shown 
in Fig. [2] where the thickness 5 of the "boundary layers" 
is to be determined as a function of Ra to satisfy Q > 0. 
With this choice of r(y) the bound will be 
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FIG. 2: Background profile with boundary layers of thickness 
< 5 < i in which r'(y) = -j-^; r {y) = for 5 < y < 1 - 5. 



Applying the horizontal Fourier transform and intro- 
ducing the shorthand D = ^ , it is evident that positivity 
of Q is equivalent to the positivity of 
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for each horizontal wavenumber k where || • || is now the 
Li^ norm on complex valued functions of ?/ £ [0, 1] and 
Re { • } indicates the real part of a complex quantity. The 
Cauchy-Schwarz and Young inequalities imply 
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Boundary conditions on Vk{y) dictate that 

Re{Dvkiy)}dy = ReK(y)}i:i;, = (21) 

so 3 j/o S (0,1) such that Re {Df)fc(?/o)} = 0. The fun- 
damental theorem of calculus followed by application of 
the Cauchy-Schwarz and Young inequalities imply 

{Rc{Dvk{y)}f =2/;^ Re{DHk{y')} Re{Dvk{y')} dy' 

< ^ {l\\Rc{DH,}f + lk^Re{Dv,}r) . (22) 



A similar pointwise bound holds for the imaginary part 
of Dvk{y) so its modulus squared satisfies 
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Thus, integrating Dv^. from to y and applying Holder's 
inequality, it is evident that 
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Because 9kiy) vanishes at y = and 1, applications 
of the fundamental theorem of calculus and Cauchy- 
Schwarz inequality yield the pointwise bounds 
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for < y < 1/2 and, for 1/2 < y < 1, 
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Restricting < 46, the task is to dominate the indefinite 
boundary layer integrals by the positive definite terms. 

The Fourier coefficients of the vertical velocity and vor- 
ticity (suppressing the time dependence) are related by 
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Integrating the modulus squared of both sides with a 
simple integration by parts implies 
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On the other hand, integration by parts and the Cauchy- 
Schwarz and Young inequalities yield 
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so that, combining (18) and (191 
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Using (|24|- 27), we conclude 
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minimize the prefactor in the bound — and minimizing 
the suitable S over k, this is satisfied by choosing S — 
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mizing wavenumber. Inserting these S and b into ( 13 1 we 
see that for Ra > 58.44 . . . (actually for Ra > ^tt^"^ 
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This exponent for the Nu-Ra upper bound scaUng, 
albeit with a prefactor 0.142, was conjectured by Otero 
from a numerical study nearly a decade ago [TB]. The 
proof here puts that result on firm analytical ground. 
The Nu-Ra and the distinguished horizontal wavenum- 
ber scaling k ^ Ra^^"* also agree with those conjectured 
by lerley, Plasting, and Kerswell following a careful com- 
bination of numerical and asymptotic analyses of the up- 
per bound problem for infinite Prandtl number Rayleigh- 
Benard convection in three spatial dimensions with free- 
slip boundaries [20] . In fact the analysis in this paper 
can be extended to that case because there is no vortex 
stretching at Pr = oo so an enstrophy balance akin to ([t]) 
is realized for free-slip boundaries |23j . 

While the rigorous bound 13 < ^ .4167 for the 
model of Rayleigh-Benard convection considered here is 
still well above that observed in most experiments and 
direct numerical simulations, it has significant ramifica- 
tions from a theoretical point of view. There are sev- 
eral theoretical predictions of Ra^^^ scaling of the heat 
transport in the "ultimate" regime of asymptotically high 
Raleigh numbers [7- 9J and the result proved here shows 
that those arguments cannot be correct without plainly 
appealing to no-slip boundary conditions or directly re- 
lying on three dimensional dynamics (or both). 

Perhaps the simplest scaling argument — making no 
mention of boundaries or boundary conditions or the spa- 
tial dimension — is the hypothesis that the physical heat 
transport is independent of the molecular transport co- 
efficients, i.e., the kinematic viscosity v and the thermal 
diffusivity k, in the fully developed turbulent regime [S]. 
This implies Nu ~ Pr^^^Ra^^^. A more physically ex- 
plicit version of the argument proceeds from the assump- 
tion that the rate-limiting process is not transferring heat 



across boundary layers into the bulk, but rather is the 
time it takes to adiabatically transport hot and cold fluid 
elements across the layer accelerated by the reduced grav- 
ity aATg neglecting frictional forces. Then the vertical 
velocity scale of rising or falling elements is ^/gaSTh and 
their heat content is ©(AT), so at sufficiently high den- 
sity of such elements the heat fiux is ~ {gahY/"^ /ST^/"^ . 
When normalized by the conductive heat fiux kAT/Zi, 
this again yields Nu Pr^^^ Ra^^^. 

More sophisticated arguments [71 19] produce the sim- 
ilar predictions. It has also been proposed that the | 
exponents will appear if the physical boundary layers are 
negligible (as might be hypothesized when Ra — > oo) 
or absent altogether. This leads to the consideration of 
"homogeneous" Rayleigh-Benard convection where the 
Boussinesq equations with a linear background profile 
are posed on a fully periodic domain. Direct numeri- 
cal simulations in three dimensions and a closure theory 
have indicated that this scaling emerges for some aspect 
ratios PH although no upper bounds on the heat 
transport can possibly exist and the genuineness of sta- 
tistical steady states is questionable for this formulation 

[251 ng. 

The Nu < Ra^/^^ bound derived here raises questions 
of precisely how the spatial dimension and the nature of 
even very thin boundary layers enter into the problem at 
high Rayleigh numbers. At least in two dimensions with 
free-slip boundaries, no matter how high the Rayleigh 
number is it is apparent that boundary layers continue 
to play a limiting role in the turbulent heat transport. 
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